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Abstract. Doppler reflectometry is considered in slab plasma model in the 
frameworks of analytical theory. The diagnostics locality is analyzed for both regimes: 
linear and nonlinear in turbulence amplitude. The toroidal antenna focusing of probing 
beam to the cut-off is proposed and discussed as a method to increase diagnostics 
spatial resolution. It is shown that even in the case of nonlinear regime of multiple 
scattering, the diagnostics can be used for an estimation (with certain accuracy) of 
plasma poloidal rotation profile. 

PACS numbers: 52.70.Gw, 52.35.Hr, 52.35.Ra 
1. Introduction 

One of widespread methods used nowadays for plasma rotation velocity measurements 
is Doppler reflectometry [H |2l Oj- This technique provides measuring fluctuations 
propagation poloidal velocity which is often shown to be dominated by plasma poloidal 
rotation velocity [3]. The method is based on plasma probing with a microwave beam 
which is tilted in respect to plasma density gradient (see figure EJ). A back-scattered 
signal with frequency differing from the probing one is registered by a nearby standing 
or the same antenna. The information on plasma poloidal rotation is obtained in 
this technique from the frequency shift of the backscattering (BS) spectrum which is 
supposed to originate from the Doppler effect due to the fluctuation rotation. 

Analytical theory of Doppler reflectometry was developed in recent papers jHEHHj, 
using analytical approach in slab plasma model, which is reliable for elongated plasma 
of large tokamaks. The linear case of probing wave single-scattering is considered there 
as well as the nonlinear process of the signal formation due to multi-scattering effect, 
which is essential for long probing ray trajectory, typical for large fusion devices. The 
diagnostics spatial and wavenumber resolution is determined and means to increase the 
method locality are discussed. 

The present paper is devoted mainly to two following topics. First of all, we modify 
the linear theory of Doppler reflectometry taking into account possible antenna focusing 
in toroidal direction, which allows the diagnostics spatial resolution to be enhanced 
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without deteriorating the poloidal wavenumber selectivity, which takes place in case 
of poloidal focusing, discussed in |3]. Secondly, we compare the diagnostics locality 
in linear and non-linear regime and dwell upon experimental evidences, allowing us to 
distinguish these two cases. 

2. Toroidal focusing in linear theory of Doppler reflectometry 

In our consideration we follow our paper taking into account possible antenna 
focusing in toroidal direction, which was not considered there. The study is made in the 
frameworks of geometrical optics (or WKB) approach, and the reader is referred to [I] 
for more accurate procedure, applied to the cut-off vicinity, where WKB approximation 
is not valid. 

We consider normalized antenna electric field in the following form 



E a (7) = e z I ^^w(x, ky, k,)f(ky, k z y k ^ 



-co 



where y, z axes denote poloidal and toroidal directions. Factor f(k y ,k z ) takes into 
account the antenna pattern describing antenna radiation in vacuum 

f(k y ,k z ) = J^ J dydzE (x = 0, y, z)<r*w-*« 

— CXD 

where Eq is vacuum antenna field, differing from E a by the absence of the reflected wave 
contribution. We consider tilted gaussian antenna pattern 



f(k y , k z ) = 2^/np y p z exp j-i 



where a possibility to provide antenna focusing in toroidal direction is taken into account. 
Corresponding parameter 1Z in case of 

^ » Pi (2) 
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has a meaning of a wavefront curvature radius at the antenna. In (|T|) K corresponds 
to the antenna tilt (JC = cj/csin??, where $ denotes tilt angle in respect of the density 
gradient). 

According to j7j radial distribution of ordinary wave electric field in WKB- 
approximation has the following form: 

W(x,k y ,k z ) = 4i 



X cos 



c^k x \Xj ky, k z j 

4 



exp 



k x (x , ky, kz) dx 



k x {x , ky, k z ) dx 



where k x {x, k y , k z ) = k 2 (x) — ky — k 2 = [oj 2 — oj 2 e (x)}/c 2 — k 2 — k 2 , the turning point 
x c (k y , k z ) is determined by the equation 

k x [x c (ky, k z ), ky, k z ] = 

and x = corresponds to the plasma border. 

The scattering signal according to reciprocity theorem |H| can be written as 



ie 



where Pi is the probing wave power and 5nn(r) is the density fluctuation with 
frequency Q. 

Using the same procedure as described in jl] we consider the turbulence to be 
slightly inhomogeneous along x direction and rotating with plasma in the poloidal 
direction, so that the density fluctuation correlation function takes the form 



(5n(x,y,t 1 )5n(x',y f ,t 2 )) = Sn 2 



x + x'\ f + °° dxdgd^ ._ , _. l2 

n(x, q,il)\ 



(2tt)* 



x exp 



ix(x — x 1 ) + iq(y — y') — if2(ti — t 2 ) 



iqv 



x + x' 



(3) 



where v (x) is the radial distribution of the plasma poloidal velocity. This allows us to 
obtain a spectral power density of the received signal in the following form [3] 



p{u s ) = (\A S \ 2 ) = P t 



dx Sn 2 (x)S(x) 



(4) 



The scattering efficiency S(x) can be shown to consist of backscattering (BS) and 
forward scattering (FS) contributions 

S(x) 



If e 2 V f + °° dg 

o 2 / —[S B s(x,q)+S FS (x,q)} 

2 \m e c 2 J 2vr 



where 



l/(-g/2,0)|- 
k 2 x (x,IC,0) 



\n [2mk x (x, /C, 0) , q, Q 



qv[x) 



m=±l 



p\ + — [Ao + mA(x)f 



-1/2 



p 4 z + — [A - K + mA(x)f 



-1/2 



(5) 
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Here m — ±1 corresponds to BS after and before the cut-off in respect of the probing 
ray propagation. The FS efficiency takes the form 



Sfs(x, q) = 2 



4 i L " 

Py + — 



A 2 



-1/2 



4 , C 
Pz + — 



x exp 



PvQM 



■1' 



where 



A(x) = - 
c 



Ao 



x c (tC,0) 



(A - n) 

2 (x,/C,0) 



-1/2 



n 



'(- 

g 2 A(x) 
2A;(x)A 



9 







g, f2 — gt>(x) 



(6) 



da;' 



A = A(0) 



k x (x',JC, 0) : 

We consider expressions (|5|) . (jBJ) from the diagnostics locality point of view. First of 
all, the locality is determined by reversed square of the radial wavenumber k~ 2 (x, /C, 0). 
This factor, corresponding to WKB-behavior of antenna electric field, underlines the 
cut-off vicinity, but for unfavorable density profiles, for example, linear or bent down 
ones, it does not suppress enough plasma periphery contribution. For the BS the second 
factor is fluctuations spectral density \n [±2k x (x, /C, 0) , q, Q] | 2 . Due to the dominance of 
long scales in the turbulence spectrum this factor underlines the cut-off vicinity, where 
k x (x, JC, 0) is small. For FS contribution the signal suppression due to the acquisition 
by antenna pattern periphery, described by the factor 



exp < — 



An 



plays an analogical role, it underlines the cut-off vicinity, where A(x) is small. 

Additional localization for BS contribution can be provided by antenna focusing to 
the cut-off, which occurs at 72. = A . If the beam is narrow enough in toroidal direction, 
so that condition (0) is satisfied for 1Z = A , the factor 

-1/2 



pi + % [A - n + mA(x)f 



-1/2 



cA(x) 



is large in the cut-off vicinity. It should be noted that, the focusing in the poloidal 
direction can not give us such an effect, due to the fact, that we can not provide 
narrow antenna beam in poloidal direction without deteriorating the diagnostics poloidal 
wavenumber selectivity, which is described by \ f (— q/2,0)\ . 

Considering toroidal focusing influence on the FS efficiency, it should be noted that 
the focusing increases the amplitude of FS signal ©, but does not improve its locality. 

The influence of the factors discussed on the spectrum of the registered signal 
is illustrated in the modelling, described in [3], and the results can be found below 
in section |U (see figure El and figure EJ) in comparison with spectra, modelled in nonlinear 
diagnostics regime. Now let us dwell upon the spectrum modelling, illustrating the 
antenna toroidal focusing influence. We consider the density profile of DIII-D tokamak 
plasma with internal transport barrier (figure Efa)) [ID]- The probing is performed 
at different frequencies and therefore with different cut-off positions. Here we take 




Figure 2. Antenna focusing influence, (a) DIII-D density profile ^Uj. (b) Turbulence 

amplitude assumed, (c) Poloidal velocity profile ( ), and velocity estimated using 

Doppler reflectometry signal frequency spectrum shift: • — using antenna focusing, 
o — without focusing. 



into account the distance between antenna and the plasma, which was assumed to be 
equal 1 m, and suppose that focusing is performed into narrow in the toroidal direction 
spot (p z ~ 1 cm) to provide condition ((2)) to be satisfied for 1Z = A . Besides that we 
take into account the turbulence suppression in the barrier region (see figure 12(b)). 

Despite the fact that density profile in the barrier region is favorable for the 
diagnostics |lj , antenna focusing makes the spectrum shift more adequate to the behavior 
of plasma velocity in the cut-off. 



3. Noninear theory of Doppler reflectometry 

In the section we review briefly the nonlinear analytical theory of Doppler reflectometry, 
which is considered in details in [HJ Ej- In case of long enough trajectory length and 
sufficient turbulence amplitude, when the following criterion 

is satisfied, where u>i is the probing frequency, 5n/n c is the turbulence amplitude, 
normalized to the density in the cut-off, x c is the distance to the cut-off, and l cx is the 
turbulence radial correlation length, we can neglect the BS during the wave propagation 
and consider multiple FS of the probing wave only. The density fluctuations in this case 
can be taken into account as a phase modulation during the probing wave propagation 
to the cut-off and backward. It can be shown that condition J7J) holds true in large 
plasma devices even at small density perturbation level 5n/n c < 10~ 2 . 
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The wave electric field is determined by Helmholtz equation 



AE + [k 2 {x] + Sk 2 {x, y, t)}E = 



where 



Sk 2 (x,y,t) 



uf6n(x,y,t) 



cr n c 

is the fluctuation of the wavenumber (here it is given for ordinary wave, but it can be 
easily written for extraordinary wave too). The electric field can be represented in the 
following form 



E(l,y) 



G[l,y\0,y o ;t]E®(yo) dy 



(8) 



where I is a coordinate along the ray trajectory and 

nu in + i / \(y 
G[l,y\0,y o ;t] = J- — -exp 



2c 



yo) 



I 



1 

n c Jo 



n[x(l'),y(°\l'),t\ dl' 



17T 

T 



Equation (jSJ) describes the transportation of initial condition from the plasma border, 
where probing antenna is situated, inside the plasma along the ray trajectory. A function 
G contains the turbulence phase shift in question, which is determined by the density 
fluctuations and represents multiple FS effect. 

According to the reciprocity theorem [HUH], the registered signal is determined as 
an integral over all plasma border of the wave reflected by the cut-off 

c 



A, 



16vr 



dyE(2A ,y)E^(y) 



(9) 



with a weight function, determined by the electric field of the acquisition 
antenna Ea\y), if we consider it as probing one. Actually, due to narrow in the 
wavenumber space antenna pattern, the main component in this signal is formed due to 
the multiple FS, which changes essentially a poloidal wavenumber of the probing wave, 
meanwhile the radial wavenumber changes the sign due to reflection off the cut-off. 

To obtain the spectrum of the registered signal and at the same time to analyze the 
diagnostics locality, as above we consider inhomogeneous turbulence, poloidally rotating 
with a plasma ©• Averaging \A S \ 2 we obtain the spectrum in question [HUH] 



S(u>) oc exp { 



21C[p- 2 + Cq 2 ) Cq 2 v 



1 2 



£(tt 2 + q 2 v 2 ] 



-i 



tq 2 v 



(10) 



Here an operator L has the meaning of the integration over distance from the plasma 
border to the cut-off, with the averaging over the turbulence spectrum 



c 2 n^ 



da; 5n (x) 



dxdgdfi 2 
(27r )3 \ n (*,q,n)\ £ 



u 2 /c 2 5{x) k~ 2 (x) 
4L, 



X c X ^ ^cxj ^ 



:n) 
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where L = [dlnn e (x)/dx| x=Xc ] _1 is the density variation scale in the cut-off. 
The integration in (fTTj) is performed with a weight function, proportional to the 
inhomogeneous turbulence amplitude, and the factor, underlying the cut-off vicinity. 

To analyze this expression for signal spectrum we consider simple case of 
homogeneous plasma poloidal rotation v(x) = v, which gives Cq 2 v = vCq 2 . In case 
of strong nonlinear regime, when antenna beam divergence is completely determined by 
the turbulence (p 2 Cq 2 ^> 1) it can be seen that spectrum frequency shift is determined 
by traditional (linear) Doppler effect 

Umax =Ui- 2KV (12) 

On the contrary, the spectrum broadening 



V 7 ^ 



is strongly influenced in nonlinear case by the turbulence amplitude and differs from 
linear one 

1/2 



lin 



dxdg df2 



\n(x,q,fl)\ fl 2 



by the factor, which can be estimated as Au/Auu n ~ 7, where 7 is determined by (J7J). 
Thus in nonlinear case the frequency spectrum width can be substantially larger than 
in linear situation, when the factor 7 is similar or less than 1. 

In case of inhomogeneous plasma poloidal rotation the spectrum frequency shift is 
actually determined by the specifically averaged rotation velocity 

VCCq 2 v 

Umax — Ui j— {16) 

p 2 + Lq 1 

In this case the frequency spectrum shift can be produced by the region with high 
amplitude of the turbulence as well as by the region with high poloidal velocity. 

The frequency spectrum broadening, which in case of homogeneous plasma poloidal 
rotation is caused by intrinsic frequency spectrum of fluctuations, is influenced here by 
additional factor associated with poloidal rotation inhomogeneity. 



Au 



t (n 2 + q 2 v 2 ) - 



tq 2 v 



Cq 1 



1/2 



4. Discussion 



At first, we consider the spectrum modelling, which was carried on according to the 
results of linear and nonlinear consideration of the Doppler reflectometry. 
In calculation the following assumptions are made: 

(i) Geometrical parameters taken correspond to Tore Supra experiments 0: u>j/c ~ 
12 cm -1 , p ~ 14 cm, $ ~ 11.5° distance to the cut-off L ~ 20 cm. 
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Figure 3. The signal spectrum evolution, (a) Assumed density profile. R denotes the 
major radius, (b) Assumed poloidal velocity profile ( -), cut-off positions (•) and 
measured poloidal velocity profile in nonlinear (•) and linear (x) diagnostics regimes, 
(c) Absolute value of frequency shift (•) and width (•) of signal spectrum related to 
different probing frequencies (plotted via corresponding cut-off position). 

(ii) We consider the same antenna for the probing and reception. The probing is 
performed at different frequencies and therefore with different cut-off positions. 

(iii) For the sake of simplicity we suppose the turbulence level to be uniform (n = 0.03) 
and its wavenumber spectra to be gaussian. The fluctuations are believed to be low- 
frequency to neglect the spectrum width in case of homogeneous poloidal velocity 
profile. 

First of all we consider plasma density profile (figure Efa)) similar to observed in 
Tore Supra [12] and step- like plasma poloidal velocity distribution (figure Efb)). The 
registered signal frequency shift (figure Etc)) is calculated using (|T3J), but, as it 
is usually done in experiment results interpretation, the measured poloidal velocity 
profile (figure Efb)) is deduced from frequency spectrum shift using equation (fT2"|) for 
traditional (linear) Doppler effect in assumption that the registered signal spectrum 
shift corresponds to velocity in the cut-off. 

It can be seen that the value of poloidal velocity measured in such a way coincides 
with assumed one in case of homogeneous poloidal velocity distribution and differs 
from it when the cut-off is situated in the region of variable velocity. In this case large 
contribution to the Cq 2 v value is made by far from the cut-off regions due to unfavorable 
bent-down density profile, leading to the obscuration of the velocity in the cut-off. Thus 
such an interpretation of the diagnostics results gives the value of the poloidal velocity 
averaged in a specific way over plasma volume. 

The frequency spectrum shift is compared in figure Efc) with spectrum width, 
calculated for £ cy ~ 2 cm. It can be seen that spectrum width can be larger then 
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Figure 4. The signal spectrum evolution, (a) Assumed density profile. R denotes the 

major radius, (b) Assumed poloidal velocity profile ( ), cut-off positions (•) and 

measured poloidal velocity profile in nonlinear (•) and linear (x) diagnostics regimes, 
(c) Absolute value of frequency shift (•) and width (•) of signal spectrum related to 
different probing frequencies (plotted via corresponding cut-off position). 

frequency spectrum shift, which is typical for Doppler reflectometry experimental 
results. Besides, figure Ofc) demonstrates the spectrum width behavior which was 
described above: the spectrum is not broadened when probing wave propagates only in 
the region with homogeneous poloidal velocity (here we neglected the broadening, which 
arises due to intrinsic frequency spectrum of the turbulence) and broadens, when the 
cut-off, which is the bound of propagation region, crosses the point, where the velocity 
changes. 

In addition, in figure Efb) we compare the diagnostics results in linear and nonlinear 
regimes. In the both cases the plasma density profile is unfavorable for the diagnostics 
spatial resolution, but one can see that nonlinear regime is more problematic from the 
results interpretation point of view. 

The second example is to emphasize the importance of plasma density profile for 
frequency spectrum formation. We consider plasma density profile (figure Ufa)) bent 
down in plasma periphery and bent up in the core. The poloidal velocity profile has high 
gradient in the 'barrier' region (figure Hfb)). It is easy to see that bent-up density profile 
underlines the cut-off contribution and improves the locality of the method. Also this 
effect is essential for linear Doppler reflectometry j3] . It leads to the fact that in case of 
favorable density profile the poloidal velocity profile measured by Doppler reflectometry 
corresponds to the certain extent to the real one. 

Let us discuss the situation when the nonlinear theory developed should be applied. 
The nonlinear regime of Doppler reflectometry diagnostics is considered in the present 
paper, when multiple FS influence is essential for registered signal spectrum formation. 
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Figure 5. Diagnostics scheme. 1 — emitting and receiving antenna, 2 — additional 
receiving antenna. 



For this situation to take place, two important criteria should be satisfied. At first, the 
turbulence amplitude is to be large enough to provide essential turbulent phase shift 
during the wave propagation, which is described by criterion J7J)- 

The second condition leading to the necessity of nonlinear theory application to 
the diagnostics results is that small-angle scattering contribution is substantial in the 
received signal. This means that small-angle scattering signal amplitude should be 
comparable or larger than BS signal, formed by linear mechanism. Multiple small-angle 
scattering contribution @ can be evaluated as 

02 r oV2p2 



(\A S \ ) ~ ^exp^ ^ \ 

4 cx c 7 [ 7 J 



The BS signal amplitude (j3J) can be estimated as 

/ 2V27T cx c 



Here |n(g)| 2 is the fluctuation poloidal wavenumber spectrum. The ratio of the signal 
amplitudes 

(\A S \ 2 ) 4 2 3/ exp{-2/C 2 4V7 2 } 



a = 



Ai™\ 2 \ PI 2 |n(-2/C)| 2 



For example, we consider turbulence spectral density \n(q)\ 2 = 4Tr£ cy [l + q 2 il y ] 3 ^ 2 , 
where £ cy ~ 2 cm, probing frequency / = ujijln = 60 GHz, tilt angle 9 = 30°, p = 2 cm. 
Then if 7 ~ 20 the small-angle scattering contribution larger than BS one: a ~ 1.5. 

To conclude, we consider the experimental evidences of linear or nonlinear regime. 
A reliable criteria, which is actually very difficult to realize in experiment, is the 
spectrum of the passed, reflected off the cut-off signal. If the probing line can be 
distinguished in the spectrum of the signal, registered by antenna 2 in figure El we deal 
with linear regime of the scattering. In other case, when the probing line can not be 
observed in the broadened spectrum, the diagnostics works in nonlinear regime. 
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Figure 6. Multi-scattering in small toroidal device. 

Another way to recognize the scattering regime is to compare BS signal spectrum 
width with one, provided by the antenna pattern due to the Doppler effect. In case of 
more broadened spectrum we deal with nonlinear regime of scattering. 

The criteria of nonlinear theory applicability discussed can be usually satisfied in 
large toroidal devices. But we should note that the nonlinear regime can be realized in 
small tokamak due to possible reflection of the probing signal off the chamber wall (this 
effect for FT-2 tokamak was investigated experimentally in [13]). It leads to the effective 
increase of the ray trajectory, which can provide the satisfaction of the criterion (|7J). 
The signal spectrum of the signal can be estimated in this case using the results of our 
consideration (flUjl. but the operator t (ITTj) should be multiplied by the quantity of the 
reflections off the plasma and wall. 

5. Conclusion 

In the present paper the Doppler reflectometry is considered in slab plasma model in 
the frameworks of analytical theory. The diagnostics locality is analyzed for linear and 
nonlinear regime. The toroidal beam focusing to the cut-off is proposed and discussed 
as a method to increase diagnostics spatial resolution. 

In nonlinear diagnostics regime frequency spectrum shift and width of registered 
backscattered signal is analyzed in dependence on plasma density profile, turbulence 
spatial distribution and spectrum and plasma poloidal velocity profile. It is 
demonstrated that the frequency shift is not influenced by turbulence absolute amplitude 
and gives an information on poloidal velocity averaged over the vicinity of the cut-off, 
the size of which depends on the density profile and turbulence distribution. 

Thus, even in the complicated situation of multi-scattering dominance Doppler 
reflectometry technique is proved to be able to give realistic information on plasma 
rotation. The spatial resolution of the diagnostics, however, suffers from transition 
to this nonlinear regime of scattering. The consideration presented allows the spatial 
resolution of the method to be analyzed for real experimental conditions, and thus the 
diagnostics results to be adequately interpreted. 
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